Almost exponential transverse spectra from power law spectra 
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Abstract We point out that exponential shape of transverse spectra can be obtained as the 
Fourier transform of the hmiting distribution of randomly positioned partons with power law spectra 
given by pQCD, which actually realize Tsallis distributions. Such spectra were used to obtain hadron 
yields by recombination in relativistic heavy-ion collisions at RHIC energies. 
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In relativistic heavy-ion collisions exponential trans- 
verse momentum distributions has been observed at low 
to mid transverse momenta, ~ 1 ^ 4 GeV Q, |^ |^ . 
This fact was often interpreted as an indication that the 
source of the particles can be described by a tempera- 
ture. This view climaxed in the "thermal model" which 
assumes a situation well described by equilibrium ther- 
modynamics 0. 

Regarding experimental data the exponential law is 
most prominent in the transverse mass spectra which fit 
the form exp((m — mT)/T) with mx = a/ m? -\- at 
mid-rapidity at SPS and RHIC energies quite well. A ten- 
dency of increasing slopes with increasing hadron masses 
in these spectra is usually interpreted as a sign of flow, 
on the basis of the relation T = Tq + m{v?) p. A flt to 
pion, kaon and proton spectra at RHIC gives Tq w 140 
MeV and (u^) « 0.5 - 0.6 at freeze out|3]- A possible 
interpretation of this fit is a common flow pattern and a 
common temperature, Tq, for all hadrons produced in the 
heavy-ion collision. This value, converted back to no-flow 
by using the above formula for pions gives about T « 210 
MeV, which is somewhat above the color deconfinement 
temperature predicted by lattice QCD calculations . 

Since there is no known hadronic process which could 
thermalize on such a short timescale at which this ther- 
mal state must have developed, it is tempting to consider 
quark level mechanisms as reason for the observed expo- 
nential spectra. 

Constituent quark and parton recombination models 
have been utilized as a picture of hadron formation in rel- 
ativistic heavy-ion collisions in recent years 

|1,E3- This 

idea roots in the constituent quark model of hadrons orig- 
inally combining mesons from two, baryons from three 
quark-like partons. In order to obtain an exponential 
hadron spectrum by this mechanism, one needs to start 
with parton spectra which are already exponential in the 
intermediate px regime of 2 — 4 GeV. At higher trans- 
verse momenta these spectra go over into power laws, as 
given by pQCD calculations. According to the respective 
regions of dominance exponential and power law parton 
distributions are simply added as a phenomenological ef- 
fort to generate the desired experimental hadron spectra. 



In this letter we present a mechanism, by which power 
law spectra of partons can be combined to give an ex- 
ponential spectrum. For this process, one needs several 
partons to recombine to a constituent quark in contrast 
to the final step of hadron recombination where only two 
or three such dressed quarks are used. In the following 
we show that the exponential can be constructed as the 
limiting distribution of cut power laws in the form 



KE) = 



1 



(1) 



with parameters a, b and c fitted to pQCD calculation 
results as in Refs.hJ. Here we modified somewhat the 
original formula by using E, the parton energy, instead 
of the transverse momentum px- At zero rapidity it is, 
however, equal to the transverse mass my, and for large 
transverse momenta to pt- 



E — co&YiyJ pPp 



(2) 



For large transverse momenta w{E) goes over into a 
power law spectrum: w ~ ab"^ \pt\~'^- 

Before turning to the combination of power laws let us 
brie fly r ecall the familiar case of the central limit theo- 
rem |l^. Short tailed probability distributions all lead 
to a Gaussian probability distribution for the sum of ran- 
dom deviates in the limit of infinitely many independent 
draws. It is easy to demonstrate this by considering uni- 
form random deviates Xi. Then w{xi) ~ 0.5 if Xi is in 
the interval [—1,1] and zero otherwise. We seek the dis- 
tribution Pn{x) of the scaled sum. 



i=l 



It is given by the n-fold integral 



(3) 



Pn{x) 



In order to derive the Gaussian in the n 
use its Fourier transform, 



Pn{.k) 
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dxe'^'' Pn{x). 
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The x-integral is easily done, and the remainder factor- 
izes into n equal contributions, 



Pn{k) = 



In the large n limit we arrive at 



lim ( 1 — — 



(6) 
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The inverse Fourier transformation results in the stan- 
dard Gaussian 

1 



Pn{x) 



(8) 
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FigH demonstrates that this limit is approached very 
quickly, already the sum of n = 3 uniform random devi- 
ates is distributed nearly Gaussian. 
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FIG. 1: Comparison of histograms of m = 200000 sums of 
n uniform random deviates in ( — 1, 1) scaled with yji/n and 
the limiting Gauss distribution exp(— a-^/2). 

For the general case one uses an arbitrary distribution, 
'w{xi) normalized to one, J w{x)dx = 1. The central 
moments of a distribution are given by higher derivatives 
of the logarithm of the Fourier transform. For those of 
the n-fold sum we have 



A:=0 



(9) 



Since log P„{k) = n\ogw{kan) when looking for the dis- 
tribution Pn{x) of the scaled sum x = a„ Xi, we arrive 
at the following scaling law for the central moments 

Cj{n) = naicj{l). (10) 

With a choice of a„ oc one assures that 

(11) 



Consequently the limiting distribution. Poo (a;), is Gaus- 
sian. 

Now we consider the dressing of constituent quarks 
with gluons. One may assume that the constituent en- 
ergy and mass coalesces from n equal, partonic amounts 
- a kind of additive parton model - but this is not nec- 
essary. The total mass in this case would he M = J2 
and each parton would carry a fraction of the trans- 
verse momentum, Xipr = 0{pT/n). In the following we 
assume that the spectral shape of the constituent par- 
tons (recombining later to hadrons) is proportional to the 
Fourier transform of the composite wave function. Fur- 
ther we assume, that such a constituent quark (e.g. an 
up quark) is composed of a bare up quark (the seed) and 
n gluons (in general, sea partons). Since these partons 
are almost massless, we consider the center of energy as 
the proper coordinate for the constituent wave function: 



(12) 



with E Xi = 1 . By equal sharing of the total energy, 
which is the case whenever the fusion happens at low 
relative momenta, it is the arithmetic mean (with 1/n 
scaled sum) of random parton coordinates. Its distribu- 
tion, Pn{x), is interpreted in the present model as the 
constituent quark wave function squared. Its Fourier 
transform is proportional to the momentum spectrum 
and is obtained as n-fold product of individual parton 
spectra, each carrying a fraction of the total energy. For 
the sake of simplicity representing all bare parton dis- 
tribution li by the same distribution, we have, in the 
general case, n contributions: 



(13) 
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with Xi — 1. Considering cut power law parton spec- 
tra as given by eq.l^ the n-fold spectrum a,t y — ra- 
pidity becomes 



XiMr 



(14) 



Since all are positive factors and their sum is normal- 
ized to one, the limiting distribution for n — > cx) enforces 
that all Xi ^ and hence results in an exponential dis- 
tribution, normalized to Poo (0) = 1 : 



Poo{pt) = exp 



f M -Mt 
\ b/c 



(15) 



lim Cj(n) — 0, for j > 2. 



This is the central result of the present letter. It 
predicts a spectral slope determined by pQCD as being 
T = b/c. From standard fits to the pQCD parton spec- 
tra one gets T w 170 — 210 MeV, depending on the flavor 
under consideration and details of the fit to the power 
law distribution from pQCD. This slope has nothing to 
do with the lattice QCD phase transition. 
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FIG. 2: exp(— |pt|c/6), the limiting distribution of cut power 
law Fourier transforms (1 + |pt|/6)~° and its finite n approx- 
imations. The value n = 20 is characteristic for the pion. 

Our result can be understood by noticing that the cut 
power law parton distribution of eq.l^ is the canonical 
Tsallis distribution (isj with a temperature T = b/c and 
g = 1 — 1/c. This distribution maximizes the Tsallis en- 
tropy with the usual canonical constraint Eiw{Ei) = 
const. The down-scaled product of n such distributions is 
still of the Tsallis form at the same temperature T = h/c 
but it belongs to the parameter q = 1 — l/(cn). In the 
n — > oo limit g — > 1 and the usual canonical distribu- 
tion is recovered with the familiar exponential spectrum. 
By using the rather large values for c usually obtained 
by fitting pQCD results, namely c = 8, one obtains the 
following series: qi ~ 0.875, q2 « 0.94, q-j, « 0.96, etc. 

It is interesting to recall here another physical applica- 
tion of the Tsallis distribution. Imperfectly thermalized 
systems, e.g. a finite number of oscillators, also lead to 
a q value not equal to one. In the case of {n + 1) oscil- 
lators, however, one arrives at a value bigger than one: 
q = \ -\- l One obtains a distribution reminiscent 

of the discretized path integral sum. 



jiTJ 



(16) 



while primordial parton distributions belong to a g value 
less than one. 



Of course, both distributions approach the Gibbs distri- 
bution as n ^ oo. 

We interpret the high power c w 8 — 9 in pQCD results 
as an accidental value simulating an "almost thermal" 
distribution with an effective q « 0.82 — 0.89 near to one. 
This explains why the observed hadron spectra at low 
and intermediate px appear nearly exponential in heavy- 
ion and many other collisions. Parton recombination to 
hadrons amplifies this effect leading to = 1 — 1/ {an) 
by n-fold recombination. The Tsallis distribution comes 
very close to the Gibbs distribution. As we shall argue 
in the following, recombination processes are faster and 
more effective with increasing phase space volume. Con- 
sequently heavy-ion collisions are the best experiments to 
study strongly interacting matter (QCD) as close to the 
canonical distribution as possible. (It is a further inter- 
esting question whether the primordial parton distribu- 
tion is described by a smaller power c in e+e^ collisions 
than in heavy-ion reactions. This fact could shed a new 
light onto thermal model interpretation of hadron pro- 
duction in e^e~ reactions 0]). A remaining question 
is that how large n has to be for P„ to look exponential. 
FigElshows this behavior with the parameters 6 = 2 GeV, 
c = 8 (causing T = 250 MeV). Here the convergence is 
slower than in the case of the sum of uniform random de- 
viates approaching a Gaussian. The value n = 20 would 
be characteristic for pion from additive mass estimate 
Af^/mg « 21.5 with mq ~ (to„ -|- md)/'2. 

Another question is in which proportions the con- 
stituent partons, used for recombination to mesons and 
baryons, may consist of contributions from clusters of 
perturbative partons of different complexity n. Our sug- 
gestion is, that one may not necessarily need to go far 
in n, i.e. to use already exponential spectra for recombi- 
nation, in order to have an effective parton distribution 
close to the usually applied one. In particular, and for 
the case of simplicity, we discuss the n = 2 case here: it 
improves already the original n = 1 distribution a lot in 
the direction of the assumed exponential. 

Fusion of a quark with a gluon and its radiation may 
be considered as dynamical processes between n = 1 and 
n = 2 parton clusters. In order to find an equilibrium 
ratio of these two type of partons one inspects a rate 
equation. The detailed balance leads to the familiar pro- 
portionality applied in recombination models: 



f2{E) = lfAE/2). 
The effective spectrum we consider is given by 



(18) 



1 



Wi 



(17) 



f{E) = Aiw{E, b, c) + A2W^{E, 26, 2c). (19) 

Inspecting the sum of thermal and original power law 
parton distribution yields used for recombination by 
Fries, et al. R,ef.|inl|. one concludes that the ratio 
A2/A1 w 0.15 gives already a surprisingly good approxi- 
mation in the relevant pT range from 6 « 1.5 GeV to 12 
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GeV (cf. FigO. We note that in the original spectra par- 
ton energy loss and a blue shift of the temperature of the 
thermal part due to radial flow have been included. For 
our figure we used the same single parton yield, scaled 
down by a factor 0.87 and added its square at pt/2 with 
weight 0.13. The corresponding Tsallis "temperature" is 
b/c = 190 MeV. 
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0.87 w1 + 0.13 w2 = Exp + w1 
wi --iX b/c = 190 MeV 
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T= 175 MeV 
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FIG. 3: The thermal and power law parton yield components 
used by Fries, et al. lToll . and the sum of the two contributions 
is compared to the 13-87% mixture of n = 1 and n = 2 parton 
clusters, using the power law distribution only. 

For realistic meson and baryon formation a further 
recombination should be applied for each constituent 
of a hadron, possibly with different flavor, in the final 
hadronization process. In a heavy-ion collision event 
many processes go on in parallel. While some quarks 
are dressing up to become constituent quarks, some 
other partons are fragmenting to lower pT hadrons di- 
rectly. The hadron pT spectra then are composed at 
least from these two components, but maybe also of "par- 
tially dressed" quarks if a moderate n in the production 
of their center of energy from random positioned gluons 
appear. These may contribute at intermediate pt to the 
hadron formation by recombination as well; their yield is, 
however, reduced compared to the ones behaving nearly 
exponentially. 

The competition between fragmentation and recombi- 
nation may in general ease the entropy problem of the 
pure recombination models, where finally one deals with 
fewer hadrons than partons. 

We analyze the entropy problem now for the final 
hadronization step. As an approximation established by 
the above discussion, we consider the recombination of 
exponential spectra from exponential precursors. The 
entropy can be obtained using the Boltzmann formula 



(the (7 = 1 case of the more general Tsallis entropy) 

S= I {~f\n,f)^^V (20) 



(27r)3 



while the number of particles are given by 

N= I /t^K (21) 



(27r)3 



The distribution f{p) we consider is proportional to 
the Fourier transform of the probabilities. Since they 
were normalized as Pn(0) = 1, here we have an un- 
known factor related to total numbers: / = aP. For 
example for exponential spectra to combine, one gets 
/""{p/n) ~ a"^^ * f{p)- This leads to the familiar coales- 
cence formula iV„ = CnN" (with C„ being proportional 
to (VT^y~"') as well. The n-cluster entropy becomes 



Sn ^ {d — nln a) Nn 



(22) 



due to (p) /T = d when considering d-dimensional mo- 
mentum exponentials. At the same time the relative 
number of n-clusters is related to the parameter a as 

The sum of all n > 1-cluster entropies is then the en- 
tropy of the recombined state, Sh, while Sq = Si is that 
of the primordial state. From the phase space integrals 
one gets 

5h = ^ 5„ = d ^ 7V„ - In a ^ (nNn), (23) 



n>l 



n>l 



n>l 



where the number of quarks redistributed in hadrons oc- 
curs as the last term in the expression: 



(24) 



n>l 



is the recombination sum rule (counting n quarks in an 
n-fold cluster). While this rule is quite strict in the re- 
combination of constituents, for dressing up quarks with 
gluons, of course, such a constraint is not valid. When- 
ever this sum rule applies, one obtains 



Sh — Sr 



^ n>l J 



(25) 



when expressing In a from the n — \ case of eq.lESl- bmce 
the recombined sum goes over n > 1 indices, and all iV„ 
arc by definition positive, ^ S'„ < S'l is unavoidable by 
the recombination of exponential spectra. With other 
words, after recombination there are fewer hadrons than 
constituents built into hadrons. The opening of several 
different flavor channels may ease this problem, but can- 
not completely solve it. The missing entropy is either 
generated by i) volume expansion (i.e. the flow must be 
generated or accelerated by the recombination process) 
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or by ii) fragmentation, when a single parton produces 
several hadrons (as it is known from e~^e~ 2-jet events), 
or by iii) the excitation of internal degrees of freedom of 
the hadrons. 

Based on these principles all hadronic transverse mo- 
mentum distributions measured at RHIC can be worked 
out as outlined in this letter. The previously assumed 
exponential parton spectra now are derived in the inter- 
mediate pt range by folding pQCD power law spectra. 
No assumption of a heat bath, a temperature or thermal 
equilibrium is required. We present this scenario as an 
alternative to, and at the same time a possible explana- 
tion for, the remarkable phcnomcnological success of the 
thermal model. 

We do not believe, however, that this result reduces 
the usefulness of heavy-ion collisions in the search for 
quark matter. Asymmetric jet quenching and the meson 
- baryon scaling of the azimuthal flow {v-z) indicate the 
presence of a collective state with properties of a highly 
absorptive, perhaps colored, medium. The present expla- 
nation of exponential transverse spectra from pQCD and 
random statistics phenomenology offers an alternative to 
the ad hoc assumption of a thermal state. In fact, the 
independence of dressing gluons when combined to a con- 
stituent parton is an ingredient of getting the exponential 
Pt distribution. This condition is probably satisfied only 
in heavy-ion collisions. Our model also calls attention to 
the fact that " statistical" is not always synonymous with 
equilibrium Gibbs thermodynamics. Certain dynamical 
phenomena, like e.g. self organizing criticality and - what 
we have just discussed - limiting parton distributions, be- 
long to a more general framework of statistical physics. 

The exponential TOT-spectra in fact approach at low 
Pt a Gaussian, exp(— p|,/(2mT)), the textbook case of 
a thermal distribution. It is interesting to note that 
several recent articles consider a Gauss distributed in- 
trinsic parton pT in the generalized parton distribution 
functions (GPD-s) ^^l- The Tsallis distribution predicts 
^{Pt) = m-T{b + vtlt) leading to {p\) « 1 GeV^, if us- 
ing a constituent mass of rriT = 0.3 GeV and the cut-off 
6 = 1.5 GeV, in good agreement with recent theoretical 
analysis of high energy pp data. 

We have seen, however, that the parameter T in this 
distribution is not related to a heat bath or another 
source of statistical fluctuations. It is combined from two 
pQCD parameters, from a cut-off b and a power c as be- 
ing T ~b/c. The "heating agent" in a general sense here 
is the effective number of independent dressing gluons 
per constituent, n, which should lead to a Tsallis distri- 
bution with q = l — \/{cn) quite close to an exponential 
at low Pt- If the transverse positions of the dressing par- 
tons are averaged as independent random variables, the 
distribution of their center of energy will have a Fourier 
spectrum, which is nearly exponential up to w 3 GeV. 
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